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We consider a metallic wire coupled to two metallic electrodes via two junctions placed nearby. A 
bias voltage applied to one of such junctions alters the electron distribution function in the wire in 
the vicinity of another junction thus modifying both its noise and the Coulomb blockade correction to 
its conductance. We evaluate such interaction corrections to both local and non-local conductances 
demonstrating non-trivial Coulomb anomalies in the system under consideration. Experiments on 
non-local electron transport with Coulomb effects can be conveniently used to test inelastic electron 
relaxation in metallic conductors at low temperatures. 



I. INTRODUCTION 

A direct relation between shot noise and Coulomb block- 
ade of electron transport in mesoscopic conductors is well 
known. In normal conductors this relation was established 
theoreticallyii^ and subsequently confirmed experimentally^. 
Later the same ideas were extended to subgap electron trans- 
port in normal-superconducting (NS) hybrids*. The latter 
results appear to provide an adequate interpretation for ex- 
perimental observations* of Coulomb effects in such systems. 

While all the above developments concern local electron 
transport and shot noise, the question arises if there also ex- 
ists any general relation between non-locally correlated shot 
noise in multi-terminal conductors and Coulomb effects on 
non-local electron transport in such systems. An impor- 
tant example is provided by three-terminal NSN structures 
which have recently received a great deal of attention in 
both experiments^^ and theory^" in connection with the 
phenomenon of crossed Andreev reflection. The latter phe- 
nomenon yields non-trivial behavior of the non-local subgap 
conductance in such structures. Further interesting features 
emerge if one takes into account electron-electron interactions. 
One can observe, for example, the sign change of the non- 
local conductance caused either by the influence of the electro- 
magnetic modes propagating along the wireii, or by positive 
cross-correlations in non-local current noise^^. Furthermore, 
positive cross-correlations in shot noise are directly linked 
to Coulomb ani-blockade of non-local electron transpor t^'^''^^ . 
Thus, a general relation between cross-correlated shot noise 
and Coulomb effects in non-local subgap electron transport 
in NSN systems turns out to be much richer than that in the 
local casei. 

In this paper we will address the impact of electron-electron 
interactions on non-local effects in normal metallic structures 
depicted in Fig. [T] Non-local properties of such systems turn 
out to be very sensitive to inelastic processes. At low temper- 
atures such processes in metallic conductors usually become 
rather weak and electrons can propagate at long distances, 
typically of order microns, without suffering any significant 
energy changes. Hence, provided voltage bias is applied to a 
mesoscopic conductor, its electron distribution function f{E) 
may substantially deviate from its equilibrium value univer- 
sally defined by the Fermi function /fI-B) = 1/(1 -f e-^/^). 
For example, low temperature distribution function f{E) 
may take the characteristic double-step form in comparatively 
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FIG. 1: (a) Schematics of the system under consideration. 
It consists of two metallic electrodes 1 and 2 coupled to a 
metallic wire of length L connecting the electrodes 3 and 4 
via the two junctions with resistances R\,R2, Fano factors 
I3\,I32 and capacitances Ci,C2. (b) Equivalent electric circuit 
of the system depicted in panel (a). 



short metallic wires attached to two big reservoirs with dif- 
ferent electrostatic potentialsi^. 

Further interesting effects emerge if one takes into account 
an interplay between non-equilibrium effects and electron- 
electron interactions. Consider, e.g., a tunnel junction be- 
tween two metallic leads. Provided the the junction resis- 
tance significantly exceeds that of the leads, the effect of 
Coulomb interaction can be modeled by introducing inter- 
actions between electrons and some linear electro-magnetic 
environmen t ^*'^^ . In this case the strength of Coulomb in- 
teraction is characterized by an effective impedance of the 
environment and the current across the tunnel junction reads 
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where fL,R{E) are the electron distribution functions in the 
left and right electrodes and P{E) is the probability to ex- 
cite a photon with energy E due to interaction between the 
junction and the environment. Provided the environment has 
a non-zero impedance and both distribution functions /l(-B) 
and are close to the Fermi function, Eq. ([1} yields 

the well known zero-bias anomaly on the I-V curve, i.e. the 
Coulomb blockade dip in the differential conductance dl/dV 
in the limit of low voltagesi^^. Furthermore, should at least 
one of the distribution functions deviate from the equilibrium 
one, the I-V curve can receive further significant modifica- 
tions. For instance, if one distribution function takes the dou- 
ble step formi^, it follows immediately from Eq. ([1} that the 
Coulomb blockade dip in the conductance should split into 
two separate dips . These dips can be - and have beeni^ - 
detected experimentally thus offering a possibility to investi- 
gate non-equilibrium effects with the aid of small capacitance 
tunnel junctions as it was demonstrated, e.g., by experimen- 
tal analysis of the impact of magnetic impurities on inelastic 
relaxation of electrons in normal metals^^.i*'. 

Despite clear advantages and simplicity of Eq. IT}, it might 
not always be convenient to employ in order to analyze com- 
bined effects of non-equilibrium and Coulomb interaction in 
metallic conductors. Indeed, the applicability of Eq. ([T} is 
restricted to junctions with very low barrier transmissions, 
i.e. the effect of higher transmissions cannot be correctly 
accounted for my means of this equation. The latter effect 
might be important, in particular if one needs to evaluate 
the non-local conductance. In addition, the function P{E) 
is usually evaluated under the assumption of thermodynamic 
equilibrium in electromagnetic environment, which effectively 
implies equilibrium electron distributions in both leads. If, 
however, the electron subsystem is driven out of equilibrium, 
self-consistent evaluation of P{E) might become a non-trivial 
problem. Furthermore, the function P{E) would in general be 
difficult to evaluate for effectively non-linear electromagnetic 
environments. 

The above complications are avoided within the kinetic 
equation analysis presented below. This approach only re- 
quires resistances of metallic leads to remain smaller than the 
quantum resistance unit h/e'^. Within the same theoretical 
framework it allows to evaluate both non-local shot noise and 
the effect of electron-electron interactions on non-local elec- 
tron transport in normal metallic conductors as well as to 
describe a non-trivial interplay between Coulomb effects and 
inelastic processes in such structures. 

The paper in organized as follows. In Sec. |lT] we outline 
our model and define the Hamiltonian of our system. In Sec. 
mil we analyze non-local correlated shot noise in the system 
under consideration. In Sec. llVl we extend this analysis taking 
into account electron-electron interactions and demonstrating 
direct relation between shot noise and interaction effects in 
non-local electron transport. A brief summary of our key 
observations is contained in Sec. [V] Some technical details 
are relegated to Appendices. In Appendix |X] we outline key 
steps of our derivation of the kinetic equation employed in 
our analysis. Necessary details of our solution of this kinetic 
equation are displayed in Appendix [Bl 



II. THE MODEL 

In this paper we will consider the system depicted in Fig. 
[1] It consists of a metallic wire of length L connected to two 
leads 1 and 2 by two small area junctions located at cc = Xj — 
(a;j,0, 0), j — 1,2 and two bulk reservoirs 3 and 4 at a; = 
and X = L [x is the coordinate along the wire). 

The system depicted in Fig. [T]is described by the Hamil- 
tonian 

H = Hi+H2 + H„i,o + HT,i+HT,2, (2) 

where 

are the Hamiltonians of the normal metals, 

+ Uix) + eV{t,x)^Xc.{x) (3) 
is the Hamiltonian of the wire and 

Ht,, =Y.I d''x[t,{x) e»^^(*' i,]^^(x)x^{x) + c.c] (4) 

are tunneling Hamiltonians describing transfer of electrons 
across the contacts with area Aj and tunneling amplitude 
tj{r). Here and below m stands for the electron mass, n is 
the chemical potential, the index a labels the spin projection, 
the potential U{x) accounts for disorder inside the wire and 
V{t, x) represents the scalar potential. The transmissions of 
the conducting channels of the junctions are related to the 
matrix elements of the tunnel amplitudes t^ ^ between the 
states belonging to the same conducting channel as follows 

t!^^ = |r(^"f = 47v\,uo /(l + Tv^i^j^o (5) 

where Vj {j = 1, 2) is the density of states in the corresponding 
terminal and is the density of states inside the wire. The 
barrier resistances _Ri and R2 and their Fano factors /3i and 
/32 are expressed in a standard way as 

i- = ¥ E 'P^-H Ti^' (1 - T^-i^') /e ■ (6) 

n n n 

A voltage bias, respectively Vi, V2, V3 and V4, can be applied 
to all four metallic terminals 1,2,3 and 4. 

In the setup of Fig. [T]one of the junctions, e.g. the junction 
2, may be viewed as an injector, which drives electron distri- 
bution function in the wire out of equilibrium. The junction 
1 may then be used as a detector for experimental investi- 
gation of nonequilibrium effects. One of the ways to observe 
such effects is to study the non-local differential conductance 
dIi/dV2 of our system. Clearly, in such kind of experiments 
the distance between the junctions should not exceed an effec- 
tive electron inelastic relaxation length Lin(r) which sets the 
scale for non-equilibrium effects in the wire at a given tem- 
perature. Thus, the setup of Fig. [T] may be used to directly 
measure Lin. 
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Finally we note that the above particular system geometry 
is chosen merely for the sake of definiteness. The key steps 
of our subsequent analysis and the results obtained from it 
remain applicable to a much broader class of systems than 
that depicted in Fig. [T] E.g., the wire may be replaced by a 
metallic lead of any shape, and ultimately all geometry spe- 
cific details can be absorbed in few elements of the conduc- 
tance matrix. 



III. CROSS-CORRELATED SHOT NOISE 

We begin with the analysis of shot noise employing the so- 
called Boltzmann-Langevin technique^"'^^ based on a kinetic 
equation for the electron distribution function f{t, E, x). Low 
frequency cross-correlated shot noise in multi-terminal metal- 
lic structures has already been studied before, see, e.g., Ref. 
[20I . Here we will briefly rederive and somewhat extend the 
corresponding results in order to illustrate the basic idea of 
the approach in a relatively simple case. In the next section 
we will extend this approach in order to include electron- 
electron interactions where more involved calculations will be 
necessary. 

The Boltzmann-Langevin kinetic equation accounts for cur- 
rent noise produced by the junctions 1 and 2 and has the form 



9/ _ 2 _ f-fF{E-eVit,x)) 

dt " Tin 



+ 



f-fF{E-ewi] 

f-fF{E-ew2] 
2e^uoR2 

rii{t,E)S{x - xi) + rj2{t,E)d{x - X2 
2evo 



■5{x — xi) 
■S{x — X2) 



'{7) 



Here D and vo are respectively the electron diffusion constant 
and the electron density of states at the Fermi energy inside 
the wire. We also introduced electrostatic potentials of the 
leads wi and 11)2 in the vicinity of the junctions 1 and 2, 



r-i, r2 < Ri,R2, 



(8) 



where the resistances of the leads fj are defined in Fig. [T]d, 
and Tin ~ D/Lf^ in the inelastic relaxation time. Note that 
here we are not going to discuss physical mechanisms domi- 
nating the process of electron energy relaxation at low temper- 
atures and simply treat rin as a phenomenological parameter. 

The potential V{t, x) should be determined self- 
consistently from the equation 



J dE[f{t,E,x) - fF{E - eV{t,x))] = 0, 



(9) 



which directly follows from the charge neutrality condition 
inside the normal metal. This charge neutrality condition in 
metals is a direct consequence of strong Coulomb interaction 
between electrons as well as between electrons and lattice ions. 
Integrating Eq. (0 over energy we obtain 

^^-DVi'^V{t,x)^ 

[wi -V{t,Xi)) {-W2 -V{t,X2)) ^, - 
+ ^ ^ o V. —0{X - Xi) + ' —d(x - X2) 



+ 



2e^voRi ' ' 2e^uoR2 

rji{t,E)S{x - xi) +rj2{t,E)S{x - X2) 



dE 



2e2 



(10) 



Note that inelastic relaxation time Tin drops out from this 
equation. 

The stochastic variables rii(t,E) and r]2{t,E) in Eqs. 
and (llOp account for low frequency fluctuations of the current 
carried by electrons with energy E through the junctions 1 
and 2 respectively. The corresponding correlators read^ 



{n,{t^,Ei)ri,{t2,E2)) 



R 



■S,jS{ti-t2)5{Ei- E2) 



X {|3,f{t^,E^,x,)[l- ME-^-ew,)] 

+ Pj [1 - .f{ti,Ei,x,)] fF{Ei - ew,) 

+ {l-P,)f{ti,E,,x,) [1 - f{t^,E^,x,)] 

+ Pj)fF{E, - ew,) [1 - fF{E^ - ewi)] ]. (11) 

Finally, no fluctuations occur at fully open contacts between 
the wire and the terminals 3 and 4. These contacts are ac- 
counted for by the boundary conditions 

f{t,E,x^Q) = fFiE-eVa), f{t,E,x^L)^fF{E-eVi). 

(12) 

Note that in the Eq. ((Tjl we have neglected the internal cur- 
rent noise generated in the wire^S. In order to justify this 
approximation, in what follows we will assume 



ri, r2, r 12, r 1,7-2 < Ri,R2, 



(13) 



i.e. we will assume the junction resistances to be much higher 
than the resistances of the metallic leads and the wire (see 
Fig. [T]d for the definition of the resistances). Thus, the task 
at hand is to solve Eqs. 0, pO|) supplemented by Eqs. (|11|) . 
(|12p and to evaluate the current noise in our system. 

As we already discussed above, the form of the distribution 
function inside the wire may essentially depend on the relation 
between its size L and the inelastic relaxation length Lin. 
Yet another relevant parameter to be compared with Lin is 
the distance between the two junctions \x2 — xi\. Provided 
inelastic relaxation is very strong, Lin ^ |2^2 — a^i| < L, the 
inelastic term in Eq. (O plays the dominant role and the 
electron distribution function / in the wire remains close to 
the Fermi function fpiE — eV{x)) with the voltage V{x) to 
be derived from Eq. (|10p . In the opposite weak relaxation 
limit L <^ Lin the inelastic collision integral in Eq. can be 
neglected. Of interest is also the intermediate limit of a long 
wire L 3> Lin but relatively weak relaxation \x2 — xi\ <^ Lin. 

We begin our analysis by defining the currents Ji and I2 
across junctions 1 and 2: 

J.(i) = ^ / dE[fF{E - ew,) - fit, E,xj)] + Si,, (14) 

where 51 j — J dE rjj{t, E) is the fluctuating current in the 
j-th junction. In the limit of full inelastic relaxation, Lin ^ 
\x2 — xi\ < L, the distribution function in the wire has the 
equilibrium form, and with the aid of Eq. (|lip we derive the 
zero frequency spectral noise power Sj — J dt{5Ij{t)5Ij{0)) , 



=Af^coth|| + (i-ft)2r, 



(15) 



where Vj = Wj — V{xj) are voltage drops across the junctions. 
Under the condition (1131) one flnds 



r2V3 + riVj 
ri+r2 



j = l,2. (16) 
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Naturally, Eq. (|15|l just coincides with the noise power for a 
perfectly voltage biased junction^^. 

Let us now consider the limit \x2 — x\\ < Lin ^ L. In this 
case, according to Eq. ([7} the electron distribution function 
f{t, E, Xj) deviates from the equilibrium form and fluctuates. 
Hence, the total current noise should acquire an additional 
contribution. In order to proceed let us establish the rela- 
tion between the distribution functions f{t,E,Xj) and the 
stochastic variables r/j. This goal can be achieved with the 
aid of the diffuson T>(t,x,x'), which is defined as a solution 
of the diffusion equation 



d_ 

dt 



+ 



1 



"Dit, X, x') — 



V{t, X, x') + S{t)5{x - x') (17) 



with boundary conditions 

V{t, 0, x) = V{t, X, 0) = V{t, L, x) = V{t, X, L) = 0. (18) 

The physical meaning of the diffuson D(t, x, x') is well known: 
It defines the probability for an electron injected into the wire 
at the point x' to reach the point x during the time t. We 
also define the Fourier transformed diffuson 

'D{lli,x,x') = j dte'^*'D{t,x,x'). 
The solution of Eq. ^ can be expressed in the form 

/(t, E,x) = j dt'd'^x'^^^^^^^^^ fp {E - eV{t', x')) 

■D{0,x,xi) 'D{0,x,X2) , , 
+ -TT^ — — JFiE - ewi) H — - — - — /f(-E - ew2) 



2e^uoRi 



+ V{t-t',x,X2)mi-t\E)]. 
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FIG. 2: Electron distribution function in the wire in the 
vicinity of the first junction, f{E,xi). The solid line shows 
f{E,xi) given by Eq. (|20p which is applicable at intermedi- 
ate values of inelastic relaxation. In this regime the distribu- 
tion function has three steps. The dashed line corresponds 
to elastic limit in which the distribution function f{E,xi) is 
defined in Eq. (|32l) . The system parameters are: T — 50 mK, 
■D{0,xi,xi)/2e^i'oRi = 0.3, 0(0, cci, a;2)/2eVoi?2 = 0.1, 
wi = 30 /iV, W2 = -50 ^iV, Va = 80 ^V, V4 = -80 ^V, 
ri/r = r2/r = 0.5. 



Here Gjj and Gni define respectively local and non-local con- 
ductances of our structure: 



J_ _ Vo{0,Xj,Xj) ^ ■Do{0,xi,X2) , . 

Rj 2e^voRi ' 2e^uoRiR2 ' ^ ' 



(19) 



where 'Do{^,x,x') is the solution of the diffusion equation 
(|17|l with Tin — >■ 00. One can equivalently write these conduc- 
tances in the from 



This general expression gets simplified in the limit e|V3 — 
V4I <^ TL/Lin and provided current fluctuations can be ne- 
glected, i.e. 771,2 — >■ 0. In this case the electric potential 
V{x) does not depend on time and slowly varies in space. 
Then one can approximately replace fp{E — eV{t',x')) by 
fp {E — eV{x)). Afterwards, employing the properties of the 
diffuson, one finds 



f{E,x) = 
0(0, X, X\ 



2e^voRi 



I5(0,a;,a;i) D(0,a;,a;2) 
2e^voRi 2e'^voR2 

'D{0,x,X2) 
2e^voR2 



■fF{E-ewi) + 



fp {E - eVix)) 



fF{E-ew2). (20) 



The non-equilibrium distribution function in this regime has 
three steps, see also Fig. [2] The first one comes from the 
distribution function of the isolated wire fF{E — eV{x)), 
while the other two steps, oc fp{E — ewj), originate from 
the junctions. Since the diffuson 'D{0,x,x') decays at dis- 
tances \x — x'\ > Lin, the distribution function acquires its 
equilibrium form far away from the junctions. 

The currents 7i and I2 can be evaluated with the aid of 
Eqs. ® and They read 



Ji = Gilt;! - Gni«2 + Sh, 
h — —Guivi + G22V2 + Sh- 



(21) 



Ri 



rir2 
ri + r2'' 



Gnl — 



rir2 



(ri + r-2)i?ii?2 



(23) 



Here we again assumed that ri2 <^ ri,r2, and ri,r2,fi,f2 <^ 
7?i,i?2. Let us emphasize that the results (|22|) . (|23|) were 
derived from Eq. pop and, hence, are not sensitive to inelastic 
relaxation at all. Besides that, general expressions (|22|) are 
not restricted to the wire geometry and remain valid for any 
shape of the leads. 

Finally, the noise terms Slj appearing in Eq. (I21|) read 



Sh = 



Sh 



dEdt' 



S{t-t') 



T>o{t - t', xi, xi] 



2e^voRi 



Vo{t-t',Xi,X2] 



dEdt' 



+ [S{t~t') 



Vo{t — t', X2,Xl 



Vi{t',E) 
(24) 



2eVoi?2 

T>o{t - t', X2,X2] 

2eH^2 



V2{t',E) 



(25) 



Clearly, they differ from the bare noise terms Sh since the 
contributions coming from electrons diffusing from one junc- 
tion to the other or returning back to the same junction are 
also taken into account. 
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We are now in position to evaluate the zero frequency noise 
power matrix 



Si' 



dt{Sh{t)SIj{0)). 



(26) 



With the aid of Eqs. and (|24I25|I we express the noise 

power for the first junction as 

^ii = RiGl, J dE{(3rf{E,Xr)[l- fp{E-ewi)] 

+ [1^ f{E,xi)]fF{E-eWi) 

+ {l-l3i)f{E,xi)[l-f{E,xi)] 

+ (1 - [1 - fF{E - ewi)] fF{E - ewi)} 



ew2)] 



+ R2GI, j dE{l32f{E,X2)[l^ friE 

+ l]2[l~f{E,X2)]fF{E-eW2) 
+ {l-l32)f{E,X2)[l-f{E,X2)] 

+ (1 - 132) [1 - fpiE - ew2)] fF(E - ew2)}. (27) 



This function is also illustrated in Fig. [2] Here we introduced 
the total resistance of the wire r — n + r2 + ri2 and assumed 
ri2 ^ r\,r2. The noise in the limit (|13|) . ri2 <^ r\,r2 and 
/3i = /32 = 1 becomes 



S\2 



— W(wi - V3) + — VK(wi - Vi) 
r r 



Gil 

+ Gr,lW(wi~W2) 
— Gnl 



— Wiwi - V3) + —Wiwi - Va) 
r r 



(33) 



(34) 



Comparing these expressions with Eqs. p9|) . pop . we observe 
that they coincide either provided V3 — V4 or in the large bias 
limit Wj — Va ^ T. Otherwise, every function W entering 
the result in the limit of strong relaxation splits up into two 
functions in the limit L <C iin. 



+ —W(w2 - Va) + —W{w2 - Vi) 
r r 



Substituting the distribution function (|20p into this expres- 
sion, assuming Gni <C Gii,G22, defining the function 



W{v) = ev coth 



2T' 



(28) 



and under the condition (|13l) . we arrive at the final result 

511 = Gii[l3iW{vi) + {l-l3i)2T] 

+ /3i Gnl Wivi - V2) + (1 - /3i) GniW{v2), (29) 

512 = -Gnl [l3iW{vi) + P2W{v2) + (2 - /3i - /32)2r](30) 



Noise power for the second junction 522 is defined by Eq. (|29p 
with interchanged indices 1 -o 2. Here we have introduced the 
effective non-local conductance 



Gnl ~ 



T>{0,x\,X2) 
2e^uoRiR2 ' 



(31) 



which, in contrast to Gni, is suppressed by inelastic relaxation. 
One has Gni ^ Gni if the distance between the junctions 
exceeds Lin, and Gni — Gni if \xi — X2\ <^ Lin- 

The first line of Eq. (|29|l just coincides with the standard 
expression for the shot noise of a mesoscopic conductor with 
the Fano factor /3i, while the second and third lines provide 
the corrections induced in the first junction by the second one. 
The origin of these corrections is simple: voltage bias applied 
to the second junction yields modifications in the electron 
distribution function in the vicinity of the first junction (cf. 
Eq. (I19|l 'l thus changing its current noise. 

Now we turn to the regime of a short wire, L ^ Lin, where 
inelastic relaxation can be fully ignored. Accordingly in Eq. 
(O we set Tin = 00 and repeat the above calculation in this 
limit. As a result, the distribution function in the wire ac- 
quires the four step shape 



fiE,x) 



Vo{0,x,xi) Vo(0,x,X2) 



2e?uaRi 



2e^uoR2 



—fF{E - eVa) + —}f{E - cVa) 
r r 



Vq{Q,x,xi) 
H TT^ — — Jf{E - ewi) 



2e^uoRi 

Vo{0,X,X2) ^ 

H 7r2 — 5 — JPyE ' 

2e-'uoR2 



ew2). 



(32) 



IV. NON-LOCAL ELECTRON TRANSPORT IN 
THE PRESENCE OF INTERACTIONS 

Until now we have ignored interaction effects and restricted 
our consideration to low frequency current fiuctuations. Be- 
low we will account for electron-electron interactions and eval- 
uate the interaction correction to the conductance matrix of 
our system. Extending the argumentsiiS^ we will demonstrate 
a close relation between Coulomb blockade of non-local elec- 
tron transport and shot noise in the system under considera- 
tion. For this purpose it will be necessary to go beyond the 
low frequency limit and allow for arbitrary (not necessarily 
slow) fluctuations of voltages Vj{t) across the junctions. In 
this regime the time and energy dependent electron distribu- 
tion function in the wire f{t, E, x) becomes ill-defined due 
to quantum mechanical uncertainty principle. This problem 
can be cured by employing the Keldysh Green function of 
electrons 



^ / N I dE 

G{U,t2,x)= I —t 



-iE(ti-t2) 



+ t2 



E,x 



(35) 



which fully describes electron dynamics at arbitrarily high 
frequencies. Applying the Fourier transformation (|35p to the 
kinetic equation ([7]) we cast it to the form^^ 

f J- + J- - I^V^ + — + i$(ti, x) - i<i>(t2, x)) G 

\Otl at2 Tin / 

I _jyg-i[*{*l>^)-*(«2,iE)] 

Tin sinh7rT(ti — 12) 
5{x-Xi_) _ire-"l'^i(*i'-'^i('2)] 



2e^voRi 

Six — X2) 
2e^voR2 
5{x — x\) 



G- 



sinh7rT(ti - 12) 

_jyg-»[</'2(tl)-02(t2)] 

sinh TvT{ti — 12) 



Vi{ti,t2)-— —V2{ti,t2). (36) 



Here the stochastic variables rij{ti,t2), which now also depend 
on two times, are correlated as follows 

{vi{ti,t2)nAt3,u)) = 
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SttT?,^ 7r2 ™o ((ti - t2) + - U) + ^2) 

_^J^^-^l<t>j(t3)-<t>j (t2)] 



_»'Tig-»['^'j ('l (U)] 

-G(t3,t2) 



sinh7rr(ti - U) 



(l-ft)(^G(tl,t4)G(t3,t2) 

_ jyg-»[<^>3 (tl (t4)l _ j J-g-ilfjij (*3)-<#>j (*2)] 

smh7rr(ii - t4) smh7rT(t3 - i2) 



In Eqs. (|36|) and (I37p we defined the fluctuating phases of 
the leads = dt'ewj[t') as well as the phase ^{t,x) = 

J^* dt' eV{t' , x), where V{t' , x) is the electric potential inside 
the wire which fluctuates both in time and in space and in- 
cludes interaction effects. 

Note that fully quantum mechanical description of inter- 
action effects in metallic conductors generally involves two 
(rather than one) quantum fluctuating phase flelds and 
"I> B (defined on the two branches of the Keldysh contour) ap- 
pearing after the standard Hubbard-Stratonovich decoupling 
of the Coulomb term in the Hamiltoniani^i^^. Provided in- 
teraction effects are sufficiently small (as is the case here, 
see below) one can effectively eliminate one of these fields, 
= <1>F — $B, and retain only the "center-of-mass" field 
$-1- = (<1>F -f $b)/2 — 5> >I>. The derivation of the kinetic equa- 
tion (I36|l in the tunnel limit /?i = /32 = 1 is presented in the 
Appendix [III Rigorous derivation of the kinetic equation (|36p 
based on the non-linear cr— model as well as its applicability 
conditions can be found in Ref. l23i . 

Now we turn to the expression for the current through the 
first junction Ii. In order to derive this expression it is nec- 
essary to solve the kinetic equation (I36|l . Technical details of 
this procedure are presented in Appendix [B] Here we directly 
proceed to the corresponding results. 

Let us first consider the limit of strong inelastic relaxation, 
L ^ Lin, and assume that the wire potential varies in space 
slowly enough, e\V3 — V4I <C TL/ Li^. In this case the current 
through the first junction acquires the form 



h 



Gil 



Vl 



Pi 

e 



dt 



e 



dt- 



sinh^ TvTt 



+ 



I3i 



sinh^ nTt 

dt'dt"Ki2(t 



Kii (t) sin[ei;it 
Kii(t) sin[e(ui - ^2)*] 



sinh^ TvTt' 



■ sin[(iii 



Gn 



V2 



P2 

e 



dt 



^i Vji', xi,X2) 
' 2e^uoRiR2 

V2)t"] 

ttT" 



sinh^ TvTt 



7^22 (t) sin[et;2t] 



+ ■ 



- J dt'dt"Kl2{t" 

sin[eii2t"]. 



^, V{t' ,Xl,X2) 

' 2e'^voRiR2 



sinh^ TiTt 

Here we have defined the response functions 



2ti 



(38) 



(39) 



which characterize the response of voltage fluctuations in the 
junction i on the current noise of the junction j. The corre- 
sponding impedance matrix Zij{u}) is defined in Appendix B, 
see Eq. HB10|I . As before, here the voltage drops v\ and V2 
are defined in Eq. p6|l . 

Repeating now the same calculation in the elastic limit L <C 
Lin, we obtain 



h 



Gil 



Vl 



Pi 

e 



dt- 



sinh-" TvTt 



Kii{t) 



(T2 Ti \ 

— sin[e('u;i - V-i)t\ + — sin[e(iui - V^)t\\ 



Gn 



V2 



Pi 

e 



dt- 



sinh^ TrTt 



K22{t) 



— sm e(ui2 
V r 



Tl 

Vs)t] + — sin[e(w2 



Vi)t]) 



Gnl 



dt- 



sinh^ TvTt 



Kii{t) sin[e(i(Ji — W2)t] 



dt' dt"K 12 {t" -t') 



TtT^ V{t',Xi,X2) 

sinh^TrTt" 2e^VQRiR2 



X sm[(TOi 
l~Pi 



W2)t"] 

dt'dt"Kl2{t" 



t') 



TtT^ V{t',Xl,X2) 

sinh^TfTi" 2e^voRiR2 



^— sin[e(ui2 — Vs)t"] + — sin[e(ui2 — Vi)t" 



(40) 



where the lead potentials wi, W2 are defined in Eq. ((H|. 

Eqs. psp . (|40|l represent the central results of this paper 
which fully determines the leading Coulomb corrections to the 
conductance matrix of our structure in both relevant limits 
of strong and weak inelastic relaxation. These results also 
allow to demonstrate a close relation between shot noise and 
interaction effects, which is now extended to include non-local 
electron transport. For example, the first line of Eq. (|38p de- 
scribes the standard - "local" - Coulomb anomaly caused by 
charging effects and related to local shot noiseii^. The next 
three lines in Eq. (|38p contain terms depending on the voltage 
difference vi — V2 and describing non-local effects. Their ori- 
gin can be traced back to the corresponding contribution to 
the shot noise in the first junction, cf. the second line in Eq. 
(|29p . Finally, the contribution in the last three lines in Eq. 
(|38p depends only on the voltage V2 and emerges from the last 
term of Eq. (|B4p oc (172). In the same way one can establish 
the correspondence between various terms in the expressions 
for the current (|40p and noise (|33p . p4p in the elastic limit. 
Perhaps we should also add that the above results remain ap- 
plicable to a much broader class of systems than that depicted 
in Fig. [1] E.g., the wire may be replaced by a metallic lead 
of any shape, and ultimately all geometry specific details can 
be absorbed in few elements of the conductance matrix. 

It is interesting to compare the results (|38l4Up with the 
predictions of the P(S)— theory ([l}. Employing the usual 
definition of the P{E) function^^ 



P^E) = J ^e'^^+JiiM^ 
Jii(t) = /" ^Re[Zii(a;)] 



(41) 

[cos Lot — 1] coth ^ + i sin uit 



and combining it with the solution of the kinetic equation ^ 
one can evaluate the current ([1} in the limit of low resistances 
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FIG. 3: Local (a) and non-local (b) difTerential conductances 
evaluated in the limit \xi — X2\ <^ Lin ^ L, Eqs. (|45p and 
(|46|l respectively. The system parameters are: T = 50 mK, 
TO = 1 ns, Rsi = Sn, Rs2 = 5n, /3i = /?2 = 1, Gn = 1 
mS, Gni = 0.1 mS. The curves at V2 = in the top panel 
and at Vi = in the bottom panel are shown in real scale, 
other curves are shifted vertically for clarity. Local differential 
conductance dli/dvi exhibits a small dip at vi = V2. Non- 
local conductance dl\/dv2 shows a much more pronounced 
peak at V2 = v\. 



of the leads h/e^rj <^ 1. Comparing the result with the Eqs. 
(|38I40|I in the tunnel limit Pi — 132 — 1, one observes that the 
P{E) approach reproduces the contributions containing the 
local response functions Kii{t), K22{t), while the corrections 
oc K\2(t) are missing. One can further verify that the latter 
corrections originate from the cross-correlation of the junction 
shot noises which are ignored in the formula (|41[) . 

In order to further specify our results it is necessary to make 
certain assumptions about the form of the kernels Kij{t). For 
typical experimental setups and at sufficiently low voltages 
and temperature it is reasonable to adopt the following ap- 
proximation for the elements of the admittance matrix of the 
environment Yij{ijj) (see Eq. (|B7P for their precise definition): 
Yii{lo) = 1/Rsi, Y22 = 1/Rs2, Y12 = Y21 = 0, where Rsi 
and Rs2 are effective shunt resistances. These resistances can 
roughly be estimated as 

Rs\ = +rir2/(ri +r2), Rs2 = ^2 + rxr2/{ri +r2). (42) 

In practice, the shunt resistances may deviate from these sim- 
ple estimates due to impedance dispersion in metallic wires 
at high frequencies^. Further assuming that Gni is small as 



compared to Vii, Y22 one finds K\2{t) ~ K2i{t) — and 
Kii{t) = e'^Rsi (1 - e-'/"") , K22{t) = e^Rs2 (l - e^'/^°) , 

where to ~ RsiCi ~ RS2C2 is the charge relaxation time 
which for simplicity is taken equal for both junctions. This 
simplification is by no means restrictive since in our final re- 
sult To appears only under the logarithm as an effective cutoff 
parameter. Under these conditions the current in the limit 
L ^ Lin (|38[) can be evaluated analytically and takes the 
form 



Ii — Gil 



vi riivi) 

egi 



~ 4nPiT 

■Gni Fl(Vl - V2), 

egi 



— Gni 



V2 



47r/32r 
eg2 



Fl{v2) 



(43) 



where we defined the dimensionless conductances of the envi- 
ronment gi = 2-K / Rsi, g2 ~ 2-k /e^Rs2 and the dimension- 
less function 



Fi{v) = Im 



/I . ev \ ^ / ^ 1 . ev 

hi * IH hi 

V27rrro 2ttT J \ 2ttTto 2ttT 



2ttT 



(44) 



Here '^(x) stands for the digamma function. Both local and 
non-local differential conductances read 



dh 

dvi 



Gil 



2/3i 
91 



dh _ „ 

T; — ^l-rnl 

OV2 



l~^Fiv2) 
92 



^Gnl—F(vi~V2), (45) 

91 

+ Gni — Fivi-V2), (46) 

91 



where we introduced another function 



F{v) = Re 
+ 



* 1-H 



1 . ev \ 

h t - "I" 

2tvTto 2-kT) 



l + *2^) 



V27rrro 2-kT J 



1 . e?; \ 

1 H h I 

2-kTto 2-kTJ 



2nT V 2%TJ 
In the elastic limit L <C Lin we find 



(47) 



li — Gil 

— Gni 



egi 
4tvP2T (r2 



( — Fi{wi - V3) + —Fi{wi - V4)) 
\ r r I 



eg2 



— Fi(w2 - V3) + —Fi(w2 ~ V4) 



— Gni t^l[w\ — W2). 

egi 

Accordingly, local and non-local differential conductances ac- 
quire the form 



dh 
dwi 



dh 

dw2 



Gil 



1 



2^1 (r2 
91 



(^j-F{wi - V3) + "^Fiwi - V4)) 



-Gni F[Vl — V2) 



-Gni 



(48) 



2P2 (r2 
92 



+ Gnl — F{wi-W2). 
91 



( — F{W2 - V3) + —F{W2 - V4)] 

\ r r / 

(49) 
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FIG. 4: Local (a) and non-local (b) differential conductances 
evaluated in the elastic limit L <^ Lin, Eqs. (|48|) and (|49p 
respectively. The system parameters are the same as in Fig. 
2. The voltage values are: V3 = -100 fj.V, V4 = 100 ^V, 
V2 — 300 /iV in panel (a) and Vi = 300 fiV in panel (b) . 



inelastic relaxation length Lin, at low temperatures and un- 
der non-zero voltage bias the electron distribution function 
acquires a characteristic double step form and the Coulomb 
dip in the differential conductance splits into two dips. This 
effect disappears provided inelastic relaxation becomes strong 
Lin < L. 

If two leads are attached to a metallic wire as it is shown 
in Fig. 1, the electron distribution function in the vicinity 
of one junction may also be driven out of equilibrium pro- 
vided electrons are injected through the second junction and 
do not relax their energies at distances shorter than the dis- 
tance between these two junctions. In this situation addi- 
tional Coulomb dip in the differential conductance appears. 

The latter configuration with two junctions also allows to 
study Coulomb blockade of non-local electron transport in the 
presence of non-equilibrium. It turns out that in this case an 
interplay between Coulomb and non-equilibrium effects yields 
more pronounced peaks in the non-local differential conduc- 
tance, see Figs. [SJa, IDs. This observation indicates that ex- 
periments on non-local electron transport in the presence of 
Coulomb effects can be conveniently used to test inelastic elec- 
tron relaxation in metallic conductors at low temperatures, as 
it has already been demonstrated, e.g., in experimentsi^'AS. 

The analysis developed here applies in the weak Coulomb 
blockade regime implying that either the resistances of metal- 
lic leads should be much smaller than the quantum resistance 
unit h/e^, or the temperature should exceed charging ener- 
gies of the barriers. In this regime there exists a transparent 
relation between shot noise and interaction effects in the elec- 
tron transporljii^. Here we extended this fundamental rela- 
tion to the non-local case, demonstrating that negative cross- 
correlations in shot noise are directly linked to Coulomb sup- 
pression of non-local conductance. This is in contrast to NSN 
structures where Coulomb anti-blockade of non-local conduc- 
tance may occur being related to positive cross-correlations 
in shot noise induced by crossed Andreev reflection. 



Local differential conductance (|45|l of a long wire with L 3> 
Lin is plotted in Fig. [3^. For a chosen set of parameters it 
is weakly affected by the second junction, although a small 
dip at vi = V2 is observed. In contrast, non-local differential 
conductance dli /dv2 is very sensitive to vi and has two peaks 
centered, respectively, at ?;2 = and V2 = «i, see Fig. ^jp. 
Fig. |4^ shows local differential conductance of a short wire, 
L <^ Lin, in which the electron distribution function does not 
relax. We observe that the conductance dli/dwi given by Eq. 
(|48|l has three dips centered respectively at wi = W2,V3,V4. 
Likewise, non-local conductance dli/dw2 defined in Eq. (I49p 
shows peaks at W2 = 101,1/3, V4, see Fig. |4j3. Comparing 
Figs. |3] and |4] we observe that the dip in dli /dvi (the peak 
in dli/dv2) occuring for strong inelastic electron relaxation 
splits into two dips (peaks) in the weak relaxation limit. 



V. SUMMARY 

Let us briefly summarize our key observations. 

We have demonstrated that Coulomb blockade corrections 
to both local and non-local conductances in metallic conduc- 
tors may change significantly provided the electron distribu- 
tion function in at least one of the leads is driven out of equi- 
librium. Provided the conductor length L is shorter than the 



Appendix A: Kinetic equation in the tunnel limit 

Let us briefly discuss the main steps of our derivation of 
the kinetic equation (|36|) . 

We start by defining the electron Keldysh Green function 

_ {x\{'ti,xi)x\{t2,X2) - x\{'t2,X2)x^{ti,xi)) 
<J — 7, l^lJ 

This function obeys the equation 



+i^{U,x)-i^{t2,X2) G 



.t*i{x)e 



■ (V'i,t (^1) x| {X2) - x\ {X2) Vii,t (^1)) 



2-KUo 

+ '2^0 ^^'^ '^'^'^ ^^'^ " '^^''^ ^^'^ ^^'^^ 

f^(j.]pi<p2(t2) ..... 



(A2) 
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Here we have defined the four-dimensional vectors Xj = 

Below we will stick to the diffusive limit in which case the 
electron distribution function remains isotropic. Then apply- 
ing the standard quasiclassical technique^*, we equalize the 
coordinates, x\ = X2 = x, and make the replacement 



2m 



+ iU{xi) -iU{x2) 



-DVi 



(A3) 



We further note that the operators ipj^f and xt in the vicin- 
ity of the barriers are not independent. They are related to 
each other via the scattering matrices of the barrier. Consider 
for simplicity the tunneling limit Tn^ ^ 1 in which case the 
corresponding transmission amplitudes in Eq. ([5]) read 

r^^) = ~27vi^U-17^ti^\ (A4) 



one can verify that it coincides with the correlator (|37|l in 
the tunneling limit /3i = 1. The pre- factors in front of the 
noise terms in the Eq. (|36p are again determined by com- 
parison with the noises of the junctions in the known non- 
interacting limit. The noise variable 772 is defined analogously. 
The operators &^ ,f]i, 7)2 may be treated as classical fluctuat- 
ing functions in the spirit of the ct— model and path integral 
formulation^. 

We finally note that the kinetic equation (|36|l can also be 
derived beyond the tunneling limit, i.e. for Tn ^ ~ 1. How- 
ever, in this general case the corresponding analysis turns 
rather complicated since the full scattering matrices of the 
barriers should be employed in Eq. (|A5[) . Without going into 
such complicated algebra here we refer the reader to Ref. [l^ 
where a general and rigorous derivation of the kinetic equa- 
tion (1361) has been carried out. 



Then we obtain 



n i = l,2 

n 

n j = l,2 



Here the superscript '° labels incoming waves unaffected by 
the barriers. In the tunneling limit considered here it suf- 
fices to identify the "incoming" operators with the full ones. 
Substituting the above expressions into Eq. (|A2|) . perform- 
ing the replacement (|A3|) and setting oc 5{x — Xj), 
after adding the phenomenological term describing inelastic 
relaxation we arrive at Eq. (|36[1 for the function G(ti, t2, x) — 
G{ti,t2, X, x) without noise terms. The pre- factors in front of 
the terms on the right hand side of Eq. 1)36^ are fixed by the 
requirement that in the absence of interactions the currents 
across the barriers have the standard Ohmic form Ij — Vj/Rj. 

The noise terms may be derived if one employs Eq. (|A2[) 
for non-averaged operator Green function 

-K _ XT x\ (fe, X2) - x\ it2, X2) Xt (^1, a^i) 

The noise operator rji is then defined as follows 

Vl{tl,t2) OC 

- i^^Ml-^ (^,,, (XO x\ {X2) - x\ (X2) ^i,t (^1)) 
^ .eri^^ (xt {Xi) {X2) - i^l^ {X2) xt (^1)) 



27r 



+ i 



etl {x)e 



-i(^>l(ti) 



2tt 

.eti(a;)e*'^i('2) 



■ (T^i.t (Xi) x| {X2) - xl {X2) ^i,t (Xi)) 



(xt (Xl) V-U (X2) - (X2)xt (Xl)). 



Evaluating the symmetrized correlator of two such operators, 

^ ivi (ii , t2)f]i {H,t4) + 7)1 (ts, ti)n-L iUM)) , 



Appendix B: Details of the solution of the kinetic 
equation 



In order to solve Eq. p6[) we make use of the same proce- 
dure as in Sec. IIIII With the aid of Eq. (|35p the expression 
for the current H14|) can be rewritten as 

^At) = -FT 1™ G{t,t,xj) 



eR-j 



+ CjVj +2TVf1j{t,t), 



sinh7vT{t - t') 



(Bl) 



where we added displacement currents recharging the capac- 
itors Cj . The solution of Eq. (|36|) takes the form 

G{t-,,t2,x) = , / dt'(fx' 



sinh7rr(ti - 12) 



■D{'-^-t',x,x') 



_jyg-»[*(ti.a!)-*(t2,iE)] 
+ — — — / dt' 



sinh7rT(ti - t2) 



V 



'tl+t2 ,1 



2e'^uoRi 



2e^voR2 



V I 



2 ^1 ? 



, ,1 , ti — t2 I ti — t2 

X 771 t H — . t - 



2 ' 2 



X »72 



t2 ^1 ti — t2 



(B2) 



where we defined (fjii) = (i>j{i) — ^{t,Xj) — J^' dt' evj{t'). 
Here we have already assumed that inelastic relaxation is 
strong, L ^ Lin, and that the the wire potential varies in 
space slowly enough, e|V3 — V4I ^ TL/Lin- 
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Combining Eqs. (|B1|I and (|B2|I we evaluate the instanta- 
neous current value in the first junction 



hit) 



dt' 



f Sjt-t') 



T>o it - t', xi,xi] 



Ri 

T>o {t - t', xi, X2) 



2e'^voRl 



V2{t'] 



+ Civi+Sh{t), (B3) 



where the noise term SIi is defined in Eq. (|24p with the 
following replacement J dErjj{t' , E) — > 2-Krjj{t' ,t'). Averaging 
the expression for the current (|B3|) . (I24|) over time we arrive 
at the following current-voltage characteristics 

h = Gll«l - GnlV2 

+ 2-KR^Gii{rii(t,t)) ^ 2TvR2Gni{v2{t,t)). (B4) 

It is important to emphasize that here the average values 
{rij{t,t)) differ from zero due to the presence of fiuctuating 
phases which account for interaction effects. 

In order to evaluate these averages it is convenient to split 
the time-dependent phases into regular and fluctuating parts, 



nit) ^ eVjt + Sipj{t), j = l,2, 



(B5) 



where the potentials Vj are defined in Eq. p6|l . In what 
follows we will assume that interaction effects remain suffi- 
ciently weak, which is the case provided either the resistances 
of metallic wires are much smaller than the quantum resis- 
tance unit, Ta ^ h/e^ , or the temperature is sufficiently high, 
T > e^/2Cj. In either case phase fiuctuations remain small, 
Sifj <^ 1, and the average {rji{t,t)) can be expressed in the 
form 

Note that fluctuating phases Sipj{t), in turn, depend on 
the stochastic variables rjj{t). In order to establish this de- 
pendence we will make use of Fourier transformed Eq. (|BH) 
which yields 

n,u = y — ClUJ —lUjGll[U})j \-tUjGn\(<^) — 

+ 2tt (— Cicj^ — iujGii{oj)^ 771, „ + 27ricjGni(cj)r;2,tj, 

e ^ ' e 

+ 27ria;Gni(i^)?7i,ij + 27r 6*21^^ — jajG22(t<j)) r\2^uj- 

Here we introduced the Fourier transform of the fluctuating 
currents ij^i^ — J dte^^^ {Ij{t) — (Ij)) and used the relation 
Svj^ui = —i(jjSipj^ui/e. The conductances Gii(tj), G22(i^) and 
Gni(i^) are again defined in Eqs. (I22|l where one should now 
substitute I?o(0, a;, a;') — >■ T>o{io,x,x'), i.e. these conduc- 
tances are expressed via Fourier transformed diffusons at a 
frequency w. From the equivalent circuit of Fig. lb we can 
also define the fluctuating currents 



J = l,2 



(B7) 



where Yij{Lu) is the admittance matrix of our structure. 
The off-diagonal elements ¥12(10) = Y2i{io) are responsible 
for cross-correlations between the junctions, which may be 
caused, e.g., by capacitive coupling between the leads 1 and 
2. Excluding the currents i^ij from the above equations we 
obtain 



dt'Kij{t~t')rj.j{t'), 



(B8) 



where the kernels Kij(t) read 

" doj e-"^' 



K.,{t) 



2ti -iuj + 

with Zij{(jj) being an effective impedance matrix 



(B9) 



/ -^^^C2-^G22(")-^y22(") 
Zii {^) = G„i (4 + ^21 (") 



A(u:) 



Gnl('^) + yi2(") \ 

(Pio) 



and 



A{co) 



{^iwCioj^ +Gii{uj) + Yii{uj)) 

X ( - ia;G2 + G22(l^) + Y22iL0)) 
-{Gnl{co) + Yl2{uj)Y. 



Combining Eqs. (|B6]) and (|B8|) . we obtain 

{Vi{t,t)) = J2 I dt'dt" Kjkit' -t") 



I Svi{t,t) „ „ 



(Bll) 



Due to causality the variable J7i(i) can only depend on the 
phases ipj(t') taken at earlier times (i.e. at t' < t), while the 
function Kij{t' —t") differs from zero only for t' > t" . Hence, 
the variable r]k{t") is independent of ipj(t'), and Eq. (|Blip 
can be rewritten in the form 



{mit,t)) = E J dt'dt"K,i{t' 



t") 



j = l,2 

5 



(B12) 



ipj—eVj t 



Here the correlator {r)\{t,t)'qi{t" ,t")) is deflned in Eq. (|37|) 
with the function G{t,t" ,x\) set by Eq. (|B2|) with omitted 
noise terms, i.e. with 771,2 = 0. The average value {r;2) is 
derived in exactly the same manner. 

Now we are in a position to evaluate the functional deriva- 
tive S {rii{t,t)rii{t",t")) /Sipj{t') from Eq. ((STJ. After a 
straightforward but rather tedious calculation one arrives at 
the result (l38l). 
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